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Eva Benckiser, 1 Maurits W. Haverkort, 1 Sebastian Briick, 2 ' 3 Eberhard Goering, 2 
Sebastian Macke, 2 Alex Frafio, 1 Xiaoping Yang, 1 ' 4 Ole K. Andersen, 1 Georg Cristiani, 1 
Hanns-Ulrich Habermeier, 1 Alexander V Boris, 1 Ioannis Zegkinoglou, 1 
Peter Wochner, 1 Heon-Jung Kim, 1 ' 5 Vladimir Hinkov, 1 '* and Bernhard Keimer 1 ' f 

1 Max-Planck-Institut fur Festkorperforschung, 
Heisenbergstrafee 1, 70569 Stuttgart, Germany 
2 Max-Planck-Institut fur Metallforschung, 
Heisenbergstr. 3, 70569 Stuttgart, Germany 
3 Experimentelle Physik 4, Physikalisches Institut, 
Am Hubland, 97074 Wiirzburg, Germany 
4 Division of Materials Science, Nanyang Technological University, 
50 Nanyang Avenue, 639798 Singapore 
5 Department of Physics, Daegu University, 
Jillyang, Gyeongsan 712-714, South Korea 

The occupation of d-orbitals controls the magnitude and anisotropy of the inter-atomic electron 
transfer in transition metal oxides and hence exerts a key influence on their chemical bonding 
and physical properties. 1 ' 2 Atomic-scale modulations of the orbital occupation at surfaces and 
interfaces are believed to be responsible for massive variations of the magnetic and transport 
properties, 3-9 but could thus far not be probed in a quantitative manner. Here we show that it 
is possible to derive quantitative, spatially resolved orbital polarization profiles from soft x-ray 
reflectivity data, without resorting to model calculations. We demonstrate that the method is sen- 
sitive enough to resolve differences of ~ 3% in the occupation of Ni e g orbitals in adjacent atomic 
layers of a LaNi0 3 -LaA10 3 superlattice, in good agreement with ab-initio electronic-structure cal- 
culations. The possibility to quantitatively correlate theory and experiment on the atomic scale 
opens up many new perspectives for orbital physics in d-electron materials. 



The electronic properties of transition metal oxides (TMOs) are determined by the interplay 
of the spin, charge, and orbital degrees of freedom of the valence electrons. Progress in 
understanding and predicting these properties relies on quantitative experimental information 
about the spatial variation of all three observables on the atomic scale. Powerful probes of 
the spin and charge densities are already available. For instance, neutron diffractometry and 
reflectometry are routinely used to determine the magnetization profiles in the bulk and near 
surfaces and interfaces, respectively. The valence electron charge density is more difficult to 
investigate, because most scattering probes couple to the total charge that is dominated by 
the core electrons. Recently, however, spectroscopic methods such as electron energy-loss 
spectroscopy 10 and soft x-ray reflectometry 11 have yielded atomically resolved profiles of the 
valence electron charge. 

The d-orbital degree of freedom is the distinguishing characteristic of TMOs compared to 
materials with valence electrons in the s- and p-electron shells. In bulk TMOs, spatial alternation 
of the d-orbital occupation ("orbital order") is known to generate "colossal" variations of 
the macroscopic properties, 1 ' 2 and the influence of "orbital reconstructions" on the physical 
properties of surfaces and interfaces is currently a subject of intense investigation. 3 " 9 At present, 
however, only the spatial average of the orbital occupation can be determined in a facile and 
quantitative manner, by means of x-ray linear dichroism (XLD), a method that relies on the 
excitation of core electrons into the valence d-orbitals by linearly polarized photons. 12 Despite 
its obvious scientific interest, quantitative experimental information about spatial variations 
of the d-orbital polarization is very limited. Methods used to determine the amplitude of 
staggered orbital order in the bulk are mostly qualitative and/or require extensive model 
calculations that add substantial uncertainties. 2 Even less information is available on orbital 
polarization profiles near surfaces and interfaces. Some evidence of orbital reconstructions near 
single TMO interfaces has been derived from comparison of XLD measurements in different 
detection modes 3 and from comparison to corresponding bulk data, 6 " 8 but these methods do 
not provide depth resolution. 

Here we focus on TMO interfaces, which are currently in the center of a large-scale research 
effort driven by prospects to control and ultimately design their electronic properties. 13 Since 
the orbital occupation determines the electronic bandwidth and the magnetic exchange interac- 
tions at and near the interfaces, detailed experimental information on the orbital polarization is 
essential for the design of suitable heterostructures. We show that quantitative, depth-resolved 
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profiles of the orbital occupation in TMO heterostructures and superlattices can be derived from 
linearly polarized x-ray reflectometry with photon energies near transition-metal L-absorption 
edges. The orbital occupation is extracted by the application of sum rules, analogous to XLD, 
so that model calculations are not required. The method is broadly applicable to TMO surfaces, 
interfaces, and multilayers and can be readily generalized to bulk diffractometry where it will 
allow quantitative measurements of staggered orbital order. 

We demonstrate the method by investigating a superlattice comprised of the paramagnetic 
metal LaNi0 3 (LNO) and the band insulator LaAl0 3 (LAO). The Ni 3+ ion in LNO has a 3d 7 
electron configuration, and the nearly cubic crystal field of the perovskite structure splits 
the atomic 3d orbital manifold into a lower-lying triply-degenerate t 2g level occupied by six 
electrons and a higher-lying doubly-degenerate e g level with a single electron. While in bulk 
LNO the two Ni e g orbitals (with x 2 -y 2 and 3z 2 -r 2 symmetry) are equally occupied, model 
calculations have suggested that the x 2 -y 2 (in-plane) orbital can be stabilized by epitaxial 
strain and confinement in a superlattice geometry (Fig. 1), and that the electronic structure 
of superlattices with fully polarized in-plane orbitals matches the one of the copper-oxide 
high-temperature superconductors. 14 " 16 This system is thus a prime candidate for "orbitally 
engineered" superconductivity. 

We have used pulsed laser deposition to grow a (4u.c.//4u.c.) x 8 LNO-LAO superlattice 
on a SrTi0 3 (STO) substrate (u.c: pseudo-cubic unit cell; see Fig. 1). Since the lattice 
constant of cubic STO is somewhat larger than the pseudo-cubic lattice constants of LNO 
and LAO (mismatch ~ 2 %), the superlattice is expected to be under tensile strain. The high 
crystalline quality and strain-state of the sample were verified by x-ray diffraction and reciprocal- 
space mapping (see the Supplementary Information). Transport 17 and optical ellipsometry 
measurements in combination with low-energy muon spin rotation 18 show that 4 u.c. thick 
LNO layer stacks in TMO superlattices are metallic and paramagnetic. The total thickness of 
247 ± 8 A of the superlattice was chosen to ensure that the sample is thin enough to clearly 
resolve total-thickness fringes, but thick enough to minimise contributions of Ni 2+ that might 
occur in the vicinity of the polar layer-substrate interface. 19 

The reflectivity measurements were performed at room temperature in specular geometry 
i.e. the momentum transfer q was parallel to the surface normal z (see the Methods section). 
The resulting data for photon energy 8047 ey far from resonance in the hard x-ray range 
(Fig. 2a) reflect the high quality of the investigated superstructure. At energies corresponding 



3 



to the Ni L 2 3 resonances at 854.7 eV and 872.2 ey we observe superlattice peaks up to the third 
order, denoted by SL (000, i = 1, 2, 3 in Fig. 2. To obtain the layer thickness and interface 
roughness of our superlattice, we fitted the q -dependent data using our newly developed 
reflectivity analysis program ReMagX 20 in a mode based on the Parratt formalism 21 (for details 
see Supplementary Information). The resulting profiles are shown in Fig. 2, and the structural 
parameters are given in the table below. The best fit is obtained when allowing the first LNO 
layer stack (buffer layer) grown on STO and the last LAO layer stack (cap layer) to be slightly 
thicker and different in roughness. We stress that our structural model is simple and gives an 
excellent description for the full set of data, measured for different energies and polarization 
geometries. 

The x-ray absorption spectra (XAS) measured in total-electron yield (TEY) and fluorescence- 
yield (FY) modes are shown in Fig. 3a. Due to the vicinity of the strong La M 4 white lines, 
the Ni L 3 line is only seen as a shoulder around 855 eV (see inset), but the Ni L 2 white line is 
clearly observed at 872 eV The FY data show a smaller intensity at the La white lines and an 
enhanced contribution at the Ni lines compared to the TEY. This is explained by self-absorption 
effects, limited energy resolution, different probing depths of both methods, and the fact that 
the topmost layer of the superlattice is a LAO layer. However, the dichroic difference spectrum 
clearly shows dips at the Ni L 3 and L 2 white line energies, which we attribute to natural linear 
dichroism (Fig. 3c). We obtained nearly identical results from FY and TEY data (not shown for 
clarity) . 

In general, natural linear dichroism reflects an anisotropy of the charge distribution around 
a particular ion. 12 To obtain a quantitative estimate of the imbalance in e g band occupation 
in our superlattice, we applied the sum rule for linear dichroism, 22 ' 23 which relates the total 
integrated intensity of the polarized spectra E \\ XiZ ) to the hole occupation n 3z 2_ r 2 and n x2 _ y2 
in the e g orbitals: 

rc 3z 2 -r 2 _ 3I E\\z 

H x 2_ y 2 4I E]]x - 1 E \\ Z 

By integrating the spectra in the range 853 -877 eV, after subtracting the La M 4 contribution 
estimated by a Lorentzian profile (Fig. 3b), we obtain n 3z 2_ r 2/n x 2_ y 2 — 1-030(5), independent 
of covalency. We define the orbital polarization 

(n 2 2 — rc 3 2 2) 

P = ) y - v ® 

\n x 2_ y 2 + n 3z 2_ r 2) 
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with n x 2_ y 2 and n 3z 2_ r 2 being the numbers of electrons and finally obtain P fa 5 %, neglecting 
covalency. In order to compare the lineshape of the polarized XAS spectra directly, we per- 
formed a cluster calculation for a Ni 3+ ion in an almost-cubic crystal field of six oxygen ions 
(Ae g =10meV) (see Fig. 3b). While the cluster calculation is expected to give good results for 
localized electrons in insulators, we expect some deviations in the calculated lineshape of the 
XAS spectra in the case of our metallic nickelate layers. Indeed, although we included reason- 
able line broadening (due to the finite temperature and experimental resolution) the measured 
lineshape is not fully reproduced, and we show the best result scaled to the experimental data 
in such a way that the integrated spectral weight is conserved. However, as pointed out above, 
the polarization-dependence of the spectra is independent of details in lineshape. To reproduce 
the measured dichroic difference, P = 5.5±2% higher x 2 -y 2 band occupation is needed, in 
good agreement with the result obtained from the sum rule. In particular, we stress that a 
higher occupation of the 3z 2 -r 2 band would result in a sign change of the XLD signal and can 
be excluded. 

While in XAS the averaged absorption of the LNO layers is measured, i.e. XAS oc 
Im (/ A iiVO +/ B iN0 ) for LNO layer stacks composed of two inner layers with scattering fac- 
tor f^ N0 and two interface layers with / B LJV0 (see Fig. 1), the intensity of the (002) su- 
perlattice reflection of a symmetric superlattice is mainly determined by the difference 
F (oo2) a- (i _ i) ^lno _ jLNoy Taking advantage of this relation, we studied the polar- 
ization dependence of the reflected intensity across the Ni L edge for momentum transfers q z 
in the vicinity of the (002) superlattice peak (Fig. 4). First, one notices that the lineshape of 
the scattering profiles is rather different compared to the XAS profiles and strongly depends 
on q z . Second, there is a pronounced polarization dependence, which only occurs at the Ni 
(L 3 and L 2 ) edge energies. In particular, the single-peak profile observed for cr-polarization 
is split into a double-peak profile for 7i-polarized photons. Although the relative intensities 
change somewhat with q z , nearly identical difference spectra are observed for both q z values. 
A qualitative comparison with the dichroic signal observed in the absorption spectra already 
shows a clear enhancement, which indicates a modulation of orbital occupancy within the LNO 
layer stack. 

In order to confirm this conclusion and extract quantitative information about the orbital 
occupation, we used the newly developed analysis tool 20 to simulate the constant-q z scans (right 
panels of Fig. 4), based on the structural parameters derived from q z -dependent reflectometry 
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and the optical constants obtained from FY-XAS (see Supplementary Information for details). To 
model the observed polarization dependence, a dielectric tensor of at least tetragonal symmetry 
has to be taken into account. Within the framework of the optical approach for multilayers, we 
implemented formulae from the magneto-optical formalism derived in Ref. 24 . We modeled our 
data with LNO layer stacks split into four unit-cell thick layers, labeled A and B in the following 
(Fig. 1). For each layer A and B we assumed a dielectric tensor with tetragonal symmetry of 
the form 



-LNO 

e tetra A(B) 



( e xx A 

V <%»J 



(3) 



with complex entries e n — e" + is" = x,z) calculated from the optical constants 5,/3 by 
— (1 — 5jj) 2 — fi 2 . and e J 2 } — 2(1 — 5jj)fijj. In the case of a homogeneous LNO stack with B and 
A layer having the same dichroism, the measured averaged dichroism was taken as input, i.e. 
e xx and e zz obtained from the XAS for E \\ x and E \ \ z, respectively. The simulated constant-q z 
scans for the homogeneous LNO stack cannot reproduce the large anisotropy observed in the 
experiment (see light blue/orange lines in the right panels of Fig. 4). We therefore considered 
the case of a modulated LNO stack with different x 2 -y 2 band occupation in layers A and B, 
keeping the averaged dichroism of 5.5±2% obtained from XLD fixed. To model the optical 
constants, we parameterized them in the following way: for the cubic case the tensor has 
diagonal entries s cuhic — ^{2e xx + e zz ). Using this and the parameter a e [0, 1] we can write 
the tensor elements: 

(C aA(B)) ;7 = (lT«)^ ± ^cubic, (4) 

with upper signs for layer A and lower ones for layer B. By varying the parameter a and 
comparing the normalized difference spectra of the constant-q z scans, we found the best 
agreement of simulations and experiment for a — 0.30±0.05 and a stacking sequence BAAB. 
This corresponds to a P B = 7±3% higher x 2 -y 2 population in the interface layers B and a 
P A = 4±1% higher x 2 -y 2 population in the inner layers A. There is a clear difference to the 
solution with stacking ABBA (not shown for clarity) . 

In order to compare our results to the predictions of ab-initio electronic-structure 
calculations, 14 " 16 we have performed LDA+L7 calculations for the particular superlattice geome- 
try studied in this work with 4u.c. LNO (for details see the Supplementary Information). Using 
equation (2), we obtain an orbital polarization of 6% (3%) for the outer (inner) Ni0 2 layers B 
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(A), in excellent agreement with the experimental result. The correspondence between the 
theoretical and experimental data also confirms that the orbital population in our system is not 
significantly influenced by structural or chemical disorder. 

We conclude that polarized resonant x-ray reflectometry in combination with the analysis 
method we have presented here is a powerful tool to accurately map out the orbital-occupation 
profile in TMO multilayers. The quantitative correspondence between theoretical predictions 
and experimental measurements we have demonstrated opens up new perspectives for the 
synthesis of TMO interfaces and superlattices with designed electronic properties, including 
but not limited to the nickelate system we have investigated here. Potential tuning parameters 
for a systematic design effort include epitaxial strain from the substrate and the constituents of 
the superlattice, the thicknesses of the individual layers, and the covalency of the chemical 
bonds across the interface. Finally, we note that the methodology introduced here can be 
readily generalized to other systems with nonuniform orbital occupation, such as surfaces and 
bulk systems with staggered orbital order, and thus has the potential to bring orbital physics in 
TMOs to a new level of quantitative accuracy. 

Methods High-quality superlattices with atomic-scale precision have been grown by pulsed-laser deposi- 
tion from LaNi0 3 and LaA10 3 stoichiometric targets using a KrF excimer laser with 2 Hz pulse rate and 
1.6J/cm 2 energy density. Both compounds were deposited in 0.5mbar oxygen atmosphere at 730°C 
and subsequently annealed in 1 bar oxygen atmosphere at 690°C for 30min. Here we investigated a 
(4u.c.//4u.c.) x8 LNO-LAO superlattice grown on an atomically flat [001] -oriented single-crystalline 
STO substrate. The crystallinity, superlattice structure and sharpness of the interfaces (roughness < 1 
u.c.) were verified by high-resolution hard x-ray diffraction along the specular rod. Reciprocal space 
mapping around the (103) STO peak shows that the strain induced by the substrate is partially relaxed 
within a total thickness of 247 A (see Supplementary Information for details). 

The resonant x-ray reflectivity and XAS measurements were carried out at UE56/2-PGM1 soft x-ray 
beam line at BESSY II in Berlin, Germany, using the advanced 3-axis UHV reflectometer described in 
Ref. 25 . The undulator beamline supplies 90% linearly cr and n polarized light (for the notation see 
sketch in Fig. 2a). Beamline settings have been chosen to result in a bandwidth of 0.57 eV at 840 eV 
(optimized intensity). The x-ray absorption have been measured with cr and n polarization of the 
incident light at an angle of incidence of 6 = 30°. While for cr polarization the measured intensity 
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corresponds directly to E \\ x, the intensity for E \\ z is deduced from I z = 4/3I n -l/3I x , with I n being 
the measured intensity with n polarization. 26 Reflected intensities have been measured with a diode 
and were corrected by the fluorescence background, measured independently with a second diode. All 
intensities have been normalized by the incoming intensity measured with a gold mesh. 

Density functional calculations were performed using the Vienna ab-initio simulation package (VASP) 
code. The projector augmented-wave 27 ' 28 method was used in the framework of the generalized gradient 
approximation (GGA). 29-32 The orbital occupation numbers were obtained by integrating the projected 
density of state within a fixed sphere on Ni x 2 -y 2 or 3z 2 -r 2 orbitals (see Supplementary Information for 
details). 
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Fig. 1: Sketch of the LNO-LAO superlattice with layer stacks of four pseudo-cubic unit cells (u.c, see 
the red box), investigated in this work. The modulation of the Ni 3d e g orbital occupation along the 
superlattice normal z is depicted by a different mixture of x 2 -y 2 and 3z 2 -r 2 orbitals and modeled with 
different scattering tensor/^ (see text). The orbital occupation imbalance has been overstated for 
clarity. 
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Fit parameters 



thickness (A) roughness (A) 



7 bilayer LAO/LNO (coupled) 15.2/15.4 3.5/2.0 
STO substrate oo 1.7 
LNO buffer layer 18.5 3.3 
LAO cap layer 18.3 3.0 

Fig. 2: Momentum-dependent x-ray reflectivity of the (4u.c.//4u.c.) x 8 LNO-LAO superlattice for 
photon energies of a, E = 8047 eV (Cu K a ; hard x-rays), b, E = 854.7 eV (Ni L 3 ), and c, E = 872.2 eV (Ni 
L 2 ) . All data have been normalized to 1 at q z = 0. The measurements in the soft x-ray range have been 
performed with a and n polarization of the incident x-rays (see sketch in a) . The solid black line is the 
best fit to the data when using the Parratt algorithm. The fitted parameters, roughness and thickness, 
are summarized in the table below. 
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Fig. 3: a, Polarization-dependent XAS spectrum (FY) across the Ni L 2j 3 edges for E \\ x (in-plane) and 
E || z (out-of-plane) polarization. In the inset we compare TEY and FY spectra for E \\ x in the full energy 
range including the La M 4 5 white lines, b, Polarized FY spectra after substraction of a Lorentzian 
profile fitted to the La M 4 line shown together with results for Ni 3+ XAS spectra with 5.5% higher x 2 -y 2 
occupation, obtained from the cluster calculation, c, difference spectra (E || x-E \\ z) calculated from the 
measured (blue points) and calculated data (orange line) shown in the middle panel. 
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Fig. 4: Energy scans of the reflectivity data with constant momentum transfer q z . Two q z values close to 
the (002) superlattice peak have been chosen from the q-dependent profiles (Fig. 2): q z = 0.409 A -1 
(top panels) and q z = 0.415A _1 (bottom panels). The left panels show the polarization-dependent 
experimental data. The right panels show the corresponding simulated curves for LNO layers with (i) 
homogeneous orbital occupation within the LNO layer stack (shifted by 4xl0~ 5 for clarity) and (ii) 
modulated orbital occupation of P B = 7±3 % higher x 2 -y 2 band occupation in the interface layers and 
P A = 4±1 % higher x 2 -y 2 band occupation in the inner layers. Experimental data and simulations are 
shown for u and n polarization of the incoming photons. To correct for the difference in absolute 
reflected intensity for a and n polarization due to the vicinity of the Brewster angle, we multiplied the 
data for n polarization with the factor cos(20)~ 2 obtained by approximation from the Fresnel formulae. 
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